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We investigate the quantum Vlasov equation with a source term describing the 
spontaneous particle creation in strong fields. The back-reaction problem is treated 
by solving this kinetic equation together with the Maxwell equation which deter- 
mines the induced time-dependent electric field in the system. The evolution of 
distribution functions for bosons and fermions is studied numerically. We found 
that the system shows a regular dynamic behavior if the back-reaction is neglected. 
But if the back-reaction is included, it is not the case and some stochastic features 
are clearly revealed in the non-equilibrium distribution function. 



1 Introduction 

In recent years much attention has beerLdavoted to the back-reaction (BR) 
problem in high-energy uarticle physics Ej~u and especially in early cosmol- 
ogy □. A flux-tube model El based on the Schwinger mechanism of the vacuum 
particle creation by a strong field is a commonly used model for the dynam- 
ical description of multiple particle phenomena. During the passed time this 
process of the spontaneous pair creation has been intensively investigal/jd for 
the case of a given external field, electromagnetic or gravitational one OB. In 
majority of the BR studies, some phenomenological source is introduced into 
a kinetic equation in a closjs analogy with the exact Schwinger result for a 
constant electric field (e.g. 111). Recaotli", a more consistent derivation of the 



source term has been done in Refs.u'LrB. In particular, an exact kinetic equa- 
tion for scalar and spinor QED with the non-Markovian source teruLfor the 
time-dependent but spatially homogeneous field was obtained Juj-Ref. a. Some 
properties of the kinetic equation of this type were studied inQ~Ell. 

A common feature of this approach is the observation of a very compli- 
cated pattern of oscillations in the density of created particles. The frequency 
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of these oscillations is of order of the zitterbewegung frequency what ev- 
idences a high particle density reached in a system. It is quite obvious that 
under such conditions the interactions of created particles should be taken into 
consideration. On the level of self-consistent mean-field, it is just the back in- 
fluence of the created particles on a formed field (the back reaction problem). 
For the thermalization process to be very actual for ultra-relativistic heavy-ion 
collisions, the account of direct particle-particle interactions is important. To 
simplify the problem in question, the relaxation time approximation Ej is used 
traditionally. 

Longer than 10-years story of studies towards these directions put more 
new questions than gave answers. This is related to high computational com- 
plexity of the equations under consideration which are integro-differential and 
highly non- linear ones. In addition, these equations need a renormalization 
procedure to remove the logarithmic divergences in the observable densities 
of energy and current. Besides, the standard relaxation time approximation 
turned out. to be too rough and should be improved for the proper treating of 
collisions Ej. 

This contribution addresses a single specific problem, namely, to study 
the influence of the BR on temporal evolution of the creation process and 
manifestation of this dynamics in observable fields and especially in particle 
distribution functions. For simplicity, we do not take into account the non- 
Abelian structure of color electric fields and concentrate on the back-reaction 
problem in application only to strong electromagnetic fields. However, the 
characteristics considered and the region of the field parameters used are of 
interest for the flux tube model applications to hadronic processes. We restrict 
ourselves to a simplest situation of a time-dependent space-homogeneous field. 

2 Basic equations 

Our approach to the BR problem is baaed on the following exact Vlasov-like 
KE for the distribution hmction f{P,t)W 

where E{t) - a strong homogeneous electric field C{P,t) is the dynamical 
source term describing the vacuum creation and annihilation processes within 



") We use the units h = c = I and the metric is chosen to be g^" = diag{l, —1, —1, — 1). 
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the Schwinger mechanism 

C{P,t) = ^(^'^'^^ r dt'WiP,t,t')[l±2fiP,t')]cose{P,t,t'). (2) 
Here use the notation the transition amphtudes 

WiPAt')^eE{t') ( -^Y ' , (3) 



uj^{P,t,t') \Pit,t') 
with the kinetic 3- vector momentum P{P±, P3), the degeneracy factor g and 

P{t,t')=P:, + e[A{t)-A{t% (4) 



c.2(P,t,i')=ei + P'(t,t'), ex = ^m2 + p2^ (5) 
0{P,t,t')^2 [ dt"uj{P,t,t"). (6) 



The KE (0) is a direct consequence of the corresponding one-particle equation 
of motion in the presence of a quasi-classical electric field = (0, 0, 0, A{t)), 
E{t) = 

For the subsequenl. calculations it is convenient to use the following local 
form of the KE 

% + ^E{t)§^^^W[l±2f]-2uu, (7) 
du ^, , du 

where was introduced two auxiliary real functionsQ u{P, t) and v{P, t) with the 
initial conditions /(to) = v{to) = u{to) = and = W{P,t,t), uj = uj{P,t,t). 

In the mean-field approximation, the distribution function f{P,t) allows 
one to find the densities of observable physical quantities. In particular, the 
conduction jcond{t) and nplarization jpoi{t) terms contribute into the electro- 
magnetic current density □ 

jin{t) = jcond{t) + jpol{t) , (8) 

(^^/(At), (9) 
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The KE (pi) should be combined with the MaxweU equation 



(11) 



which closes the set of equations for the BR problem. We assume that a 
particle-antiparticle plasma was initially formed due to some external field 
Eex (t) excited by an external current jex (t) ■ The internal field and current are 
noted as Ein{t) and jin{t). So, we have 



E{t) = Emit) + Eex{t) , 3tot{t) = J^n{t) + jex{t) . 



(12) 



It is well known that vacuum expectation values of type (^) can have 
ultra-violet divergences and neecLsome regularization procedure. We use here 
the method suggested in paper EJ. To regularize different observables (cur- 
rents, energy density etc.) it is necessary to fulfill some subtractions of the 
relevant counterterms from the every regularized function /, v and u. These 
subtraction terms are constructed as coefficients of the asymptotic expansion 
of corresponding functions in series over the power of uj~^{P). The leading 
terms of such expansions can be easily found from Eqs. (^ 
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^ '^UJ^ \P3 
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The conduction current (^) is regular one while the polarization current (g|) 
contains the logarithmic divergence. For its regularization it is enough to 
fulfill one subtraction w ^ i> — ti^ in (H) that can be interpreted as the charge 
renormalization. As the final result, the regularized Maxwell equation can be 
written in the following form (it is implied here that th|C-|Coupling constant e 
and the fields Em and Eex have been renormalized too)ll2l: 



Ein — 



47^3 



9-1 



8w4 



2(3-1)' 



(14) 



This regularization procedure is mostly adequate to the method based 
on reducing the KE to the system of partial differential equations (0). This 
approach leads to the numerical results (seeEJ) which are quite^cpnsistent with 
the results obtained with the adiabatic regularization schemeHfl 



3 Numerical results 

The distribution function of a strongly non-equilibrium state of particle-antiparticle 
plasma is investigated numerically for the cases with {Em{t) ^ 0) and with- 
out {Em{t) = 0) taking into account the BR mechanism. Various shapes of 
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the external field impulse were studied. It turned out that the system behavior 
is weakly sensitive to a particular shape of the impulse and below we present 
here the calculation results only for the impulse of the Narojny-type 



AeAt) - ^0 b [tanh(t/&) + 1], E^^) = Aq cosh-^ (t/b) . (15) 



The parameters of this potential are chosen in accordance with conditions of 
the flux-tube modelcl. In particular, the coupling constant is taken as a rather 
large value, = 4 and the used dimensionless variables are : t tm, P 
P/m,A —>■ A/m. The initial impulse is characterized by the width bm — 0.5 
and amplitude Aq/tti = 7.0. 

The distribution functions of bosons and fermions are shown in Fig. |l| ne- 
glecting the BR mechanism. Because in this case their momentum dependence 
is determined only by the transition amplitudes (^, these distributions are 
smooth. Such behavior is regarded as a regular one. A valley region in the 
boson distribution function arising near small values of P3 is caused by the 
linear Ps-dependence of the amphtude (||). 



Figure 1: Momentum distribution of produced bosons (left)and fermions (right) at time 
t = 0.05. The influence of the back reaction is not included. 

When the BR mechanism is taken into account, the regular momentum 
dependence of the distribution function is destroyed (Figs. |[^; the result like 
that in Fig. | was obtained previously inn but in the framework of different 
approach) . 

At the same time Fig. ^ and Fig. ^ demonstrate the existence of periodic 
temporal behavior of the distribution function. Two-dimensional representa- 
tion in Fig. ^ clearly shows how a 'dog-brush' structure of the distribution 
function along the P3 axis is combined with the periodic time structure along 
the time axis. 
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4 Conclusion 

Thus, the BR equations generate some large-scale structure on the background 
of small-scale multi-mode complex dynamics. The small-scale trembling is a 
manifestation of vacuum oscillations. Trembling frequency is increasing with 
time. The smoothed initial distribution function in Fig. ^ corresponds to the 
external field impulse. Large-scale wiggles in the distribution function is a 
consequence of self-organization of the system Ej due to the growth of the 
collective plasma oscillations. 
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The presented results provide some evidence that the inclusion of the BR 

mechanism into consideration of the pair creation gives rise to stochastic be- 
havior of the system. It is not excluded that a dynamical chaos can be found 
also at other values of the parameter eE, however this needs additional nu- 
merical investigations. It is quite possible that the revealed stochastic features 
of the vacuum pair creation is a source of the statistical behavior observed in 
multiple particle production. 
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